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ABSTRACT 

Starting out from the most general, gauge invariant and renormalizable scalar potential of the 
R-parity violating MSSM and performing a calculable rotation to the scalar fields we arrive at 
a basis where the sneutrino VEVs are zero. The advantage of our rotation is that, in addition, 
we obtain diagonal soft supersymmetry breaking sneutrino masses and all potential parameters 
and VEVs real, proving that the MSSM scalar potential does not exhibit spontaneous or explicit 
CP-violation at tree level. The model has five CP-even and four CP-odd physical neutral scalars, 
with at least one CP-even scalar lighter than Mz- We parametrise the neutral scalar sector in 
a way that resembles the parametrisation of the R-parity conserving MSSM, analyze its mass 
spectrum, the coupling to the gauge sector and the stability of the potential. 
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1 Introduction 



That none of the terms in the Standard Model (SM) violate lepton number (L) is not 
due to an imposed symmetry, but merely reflects the fact that all such combinations of 
SM fields are ruled out by consideration of gauge invariance and renormalisability [1] . For 
supersymmetric extensions of the SM this is no longer true. In the Minimal Supersym- 
metric Standard Model (MSSM) [2], lepton number violating terms (and baryon number 
(B) violating terms) appear naturally, giving rise to tree-level processes, proton decay for 
example, which are already strongly constrained by experiment. Either, bounds can be 
set on Lagrangian parameters, or a further discrete symmetry can be imposed on the La- 
grangian, such that these processes are absent. The discrete symmetry most commonly 
imposed is known as R-parity {Rp) [3,4]. Under R-parity the particles of the Standard 
Model including the scalar Higgs fields are even, while all their superpartners are odd. 
Imposing this symmetry has a number of effects. Firstly, any interaction terms which vio- 
late lepton number or baryon number will not appear. Secondly, the decay of the lightest 
supersymmetric particle (LSP) into SM particles would violate Rp ; the LSP is therefore 
stable. The sneutrino vacuum expectation values (VEVs) are zero; without extending the 
MSSM field content, spontaneous generation of Rp violating terms is phenomenologically 
discounted [5]. 

If Rp conservation is not imposed, fields with different Rp mix [6-12]. In particular, 
the neutrinos will mix with the neutralinos and the sneutrinos will mix with the neutral 
scalar Higgs fields; all five complex neutral scalar fields can acquire vacuum expectation 
values. Minimising this ten parameter potential in general is not straightforward, it is more 
convenient to simplify the system by choosing an appropriate basis in the neutral scalar 
sector. As one of the Higgs doublets carries the same quantum numbers as the lepton 
doublets (apart from the non-conserved lepton number), it is convenient to introduce the 
notation Ca — (Hi, Li) where Hi and Li are the chiral superfields containing one Higgs 
doublet and the leptons, respectively {a — 0, ... ,3 and i — 1, ... ,3). Furthermore, starting 
from the interaction basis, we are free to rotate the fields and choose the direction corre- 
sponding to that of the "Higgs" field. Assuming that the system defining the five complex 
vacuum expectation values of the fields was solved, four complex VEVs Va would define 
a direction in the four dimensional {Hi, Li) space. One can then choose the basis vector 
which defines the Higgs fields to point in the direction defined by the vacuum expectation 
values. We refer to this basis, in which the "sneutrino" (as we call the fields perpendicular 
to the "Higgs" field) VEVs are zero, as the vanishing sneutrino VEV basis [13-15]. This 
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basis has the virtue of simphfying the mass matrices and vertices of the theory and thus is 
better suited for calculations. 

Basis independent parameterisations can be chosen which explicitly show the amount of 
physical lepton number violation [16-18]. Values for physical observables such as sneutrino 
masses and mass splitting between CP-even and CP-odd sneutrinos have been derived in 
the literature in terms of these combinations but usually under some approximations (for 
example the number of generations or CP-conservation). We find this procedure in general 
complicated for practical applications and we shall not adopt it here. 

Instead, we present in the next section a calculable procedure for framing the most 
general MSSM scalar potential in the vanishing sneutrino VEV basis. An advantage of 
our procedure is to obtain a diagonal "slepton" mass matrix, two real non-zero vacuum 
expectation values and real parameters of the neutral scalar potential in the rotated basis. 
The latter proves that the neutral scalar sector of the most general R-parity violating MSSM 
exhibits neither spontaneous nor explicit CP- violation in agreement with [19]. In Section 
3, the tree- level masses and mixing of the neutral scalar sector is investigated. Using the 
Courant-Fischer theorem for the interlaced eigenvalues, we prove that there is always at 
least one neutral scalar which is lighter than the Z-gauge boson. We present approximate 
formulae which relate the Higgs masses, mixing and Higgs-gauge boson vertices of the R- 
parity conserving (RPC) case with the R-parity violating (RPV) one. In Section 4, the 
positiveness of the scalar mass matrices and stability of the vacuum is discussed. 



2 Basis choice in the neutral scalar sector 

In this section we develop a procedure connecting a general neutral scalar basis with the 
vanishing sneutrino VEV basis, the latter being more convenient for practical applications. 
The most general, renormalizable, gauge invariant superpotential that contains the minimal 
content of fields, is given by 



>V = eat 



+ 7;^^y^K,kU^DyDl, (2.1) 



where Q^^, -Df , f/f , Ei, H^, are the chiral superfield particle content, z = 1, 2, 3 is 
a generation index, x = 1,2,3 and a = 1, 2 are SU (3) and SU (2) gauge indices, respectively. 
The simple form of ()2.H1 results when combining the chiral doublet superfields with common 
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hypercharge Y = — | into £"^q 3 = (iff, -f/^^i 2,3)- f^a is the generalized dimensionful /i- 
parameter, and \af3k, ^'ajk^ Kjk^ (Xu)ij are Yukawa matrices with eab and e^yz being the 
totally anti-symmetric tensors, with eu = 6123 = +1- Then the five neutral scalar fields, 
i'La,h2 from the SU{2) doublets, C a — {yLai ^La)'^ and H2 — (^2"? ^2)"^' form the most 
general neutral scalar potential of the MSSM, 
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(2.2) 



where general complex parameters 5q, an hermitian matrix ('^^)^^ and m^^^ all arise from 
the supersymmetry breaking sector of the theory. The last term in ()2.2|1 originates from 
the D-term contributions of the superfields H2- Defining 



one can rewrite the potential in ()2.2|) in a compact form as 

Ke.t.ai = ^>L^>L/3 + ml hi* hi - (b^hahl + H.c) 

+ l{9' + 9l)Khl-ul^9,j' . 



[2.3] 



[2A) 



In order to go to the vanishing sneutrino VEV basis, we redefine the "Higgs-sneutrino" 
fields 



La — ^aP^Lp 1 



where U is a 4 x 4 unitary matrix 



U = Vdiag(e^'^-) Z 



:2.5) 



:2.6) 



being composed of three other matrices which we define below, V unitary and Z real 
orthogonal. The potential in the primed basis becomes. 



M'r Z 



a/3 



I 2 7,0* 1,0 



{h'Z)J,X + H.c] + + gl) (hl*hl - u'lAa)' , (2.7) 



where 



diag(e-*'^'^) (A<y Vdiag(e^'^") 



V diag(e 



i4>a 



(2. 



The unitary matrix, V, is chosen such that yM.'cj is real and diagonal - the hat (") is used 
to denote a diagonal matrix. The phases (pa are chosen such that b'^ is real and positive 
[they are equal to the phases of (b'^V)^]. The minimisation conditions for the scalar fields 
are now derived, to obtain conditions for the vacuum expectation values. 



dV 



dp 



La 

dV 



M' 



V 



L/3 



0* 



1 

4^ 



2 --X9^^9l)Wh\ 



m 



lhl-{b'Z)J*^ + \{g' + gl){hl*hl 



where "vac" indicates that the fields have to be replaced by their VEVs, 



in.) = ^ 







:2.9) 



(2.10) 



The U{1)y symmetry of the unbroken Lagrangian was used to set the phase of Vu to zero, 
however, at this stage all other vacuum expectation values will be treated as complex 
variables. By combining Eqs. ()2.9|2.10|) we obtain 



' M' 



(2.11) 



mlvu - ib'Z)avl + -(fi-^ + gj)ivl - V*V^) Vu 







(2.12) 



In a general basis, it is difficult to solve the above system with respect to the VEVs without 
making some approximations, for example assuming small "sneutrino" VEVs [11]. In order 
to simplify calculations we would like to find a basis where the "sneutrino" VEVs vanish, 
Vi = V2 = Vs = 0. In other words, we are seeking an orthogonal matrix Z, such that the 
following equation. 
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(2.13) 



holds. If the above system is satisfied, then a solution with zero "sneutrino" VEVs exists. 
The other solutions, with non-vanishing "sneutrino" VEVs will be discussed later. In 
Eq. dni, 

1 



4 



(2.14) 



is the Z-gauge boson mass squared. It is obvious that when vi = 0, Vq is real. It is now 
useful to define 



tan/5 = — . 



(2.15) 



To determine Z, multiplying ()2.13|) by Z^^-, summing over a and solving for Z^o, yields 

6^tan/3 



For given set of model parameters, Z^o depends only on tan/3 which we can now fix by 
solving the orthonormality condition, 

'^ZaoZaO = ^ ^ -2 = 1- (2-17) 



a=0 a- 



This equation can be easily be solved numerically for any given set of model parameters. 

It is worth noting that when 6j = and using notation more typical for this case, 
6q = m^2) = "^ii Eq- ()2.17j) reduces to one of the standard RFC MSSM equations 

for the Higgs VEVs: 



ml^Vd = Vu 



(2.18) 



For some parameter choices Eq. ()2.17|) may admit multiple solutions for tan/3. Each 
of the possible tan/3 specify a different basis, and each of these bases has one solution of 
the minimisation conditions with vanishing "sneutrino" VEVs. The subtlety highlighted 
earlier is the following: all possible solutions of the minimisation conditions can be found 
in each basis, so, in general, each basis contains a number of extrema equal to the number 
of possible solutions for tan/3. Hence, a solution with f « = in one basis, is a solution with 
fj 7^ in another basis. The important point to note is that by considering all possible 
values of tan/3, and selecting the value which corresponds to the deepest minima for the 
solution with vanishing sneutrino VEVs, all the solutions will have been accounted for, and 
the vanishing sneutrino VEV basis will have been determined correctly. The value of the 
potential at the vacuum, in terms of tan (3 is given by 

Mi /tan2/3-l\2 

The obvious conclusion from the equation above is that the deepest minimum of the po- 
tential is given by the solution for tan/3 or cot (3 which is greatest. 

Knowing tan (3, one should fix rn^ using Eqs. ()2.12f2. 14112. lf)|l (again in the analogy with 
RFC MSSM where m| is usually given in terms of M^, tan j3). Namely 

™l^Z.o*icot^-iM|!5g^. (2.20) 
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In this way is chosen to give the correct value of the the Z-boson mass. 

Only the first column of the Z matrix, Zao, is defined by Eq. ()2.16p . The remaining 
elements of Z must still be determined. Having fixed the first column of the matrix, the 
other three columns can be chosen to be orthogonal to the first column and to each other. 
This leaves us with an 0(3) invariant subspace, such that the matrix Z is given by 



O 




(2.21) 



where 



o 



00 

Z30 



/ zIq + z|o + z'Iq 










ZqoZk) 






a/ ^20 +-^30 




Z30 


•Zoo -^20 


^30 




V^10+^20 + -^30 
^10-^20 


V^l + ^20 + 

ZooZ^o 


^30 


V^lo 


^30 V ^1 + ^20 + ^30 
^10-^30 


V-^2 0+^30 
Z20 


a/^io+^Io'^ 


^30 




^3 V ^1 + ^2 + ^3 





(2.22) 



and X is an, as yet, undetermined 3x3 orthogonal matrix determined by three angles. This 
remaining freedom can be used to diagonalise Z^ (j^'c 

2 ~ 



i.e. the (real symmetric) 



"sneutrino" part of the Z'^ y^'cj ^ matrix, with entries (M~)j. We have now accomplished 
our aim of finding the matrices V and Z which, after inserting into potential of Eq. ()2.7|1 
and dropping the primes, reduce the scalar potential to the form 



V , 

* ncutr; 



Bai>Lahl + H.C 



+ li9' + 9l){hThl 



where 



(M^ 



Z^{M''c]Z 



and Bn 



af3 



(2.23) 



(2.24) 



with yAd' and h' given by Eq. ()2.8|) . In this basis the matrix M~ adopts a particularly 
simple form 



Bo tan P - iM| cos 2p Bj tan /? 
Bi tan/5 



(2.25) 



where there is no sum over i in the down-right part of the matrix. Notice that we did not 
only succeed to self consistently go to a basis where the sneutrino VEVs are zero, but also 
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we managed to have the sneutrino masses (M~)j diagonal and all the parameters of the 
scalar potential in Eq. ()2.23|) real. 

As a byproduct of our procedure, we denote here that the potential of Eq. ()2.23p exhibits 
neither spontaneous nor explicit CP-violation at the tree level. The latter is in agreement 
with the results of Ref. [19] following a different method. Of course, the parameters fia of 
the superpotential and the soft supersymmetry breaking couplings stay in general complex. 
The result that the neutral scalar potential is CP invariant can also be seen directly from 
Eq. ()2.4jl . By forming the complex basis {i^Lay^2*) ^^le first line of the potential can be 
rewritten as a matrix; a rotation can then be performed such that the matrix is real 
and diagonal. After the rotation, the second line, being the contribution from D-terms, 
contains complex parameters in general, but the rotation matrix can be chosen such that 
these phases are set to zero. 

A question arises when we include high order corrections to the potential. Then the 
vanishing "sneutrino" VEVs will be shifted to non-zero values by tadpoles originating, for 
example, from the CQD contribution in the superpotential ()2.1|) . The "sneutrino" VEVs 
maybe set back to zero by a renormalization condition such that a counterterm for these 
VEVs set their one particle irreducible (IPI) tadpole corrections to zero. 

To conclude, it is worth making a remark about the sign of Bq. As is clear from the 
form of Eqs. ^2MMM^ if (M'I) - | M| > for all a, Bo is always positive 

in the vanishing sneutrino VEV basis. 

3 Parametrising the neutral scalar mass matrices 

The neutral scalar sector of the R-parity violating MSSM is in general very complicated. 
This is due to the fact that the scalars mix through the lepton number violating terms 
proportional to Bi, Ml and unless all of these parameters and VEVs are real one has 
a 10 X 10 matrix to consider. However, for any given set of model parameters, one can 
always perform the basis change described in the previous section and arrive to the potential 
defined by Eq. ()2.23j) . with only real parameters. Consequently, the physical neutral scalars 
are, at the tree level, exact CP-eigenstates. This implies that the neutral scalar mass 
matrix decouples into two 5x5 matrices, one for the CP-odd particles and one for CP- 
even. In the same manner as in the R-parity conserving MSSM, once quantum corrections 
are considered, the CP invariance will generically be broken [20]. 

Ultimately, one would like to parametrise the scalar sector resulting from the potential 
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in ()2.23|) with as few parameters as possible in order to make contact with phenomenology. 
These parameters in the case of the R-parity conserving MSSM are: the physical mass of 
the CP-odd Higgs boson 

2 So 



Ml 



sin 2/? 



(3.1) 



and tan/5. An advantage of the form of potential in Eq. ()2. 2312. 2412. 25|) is that, Ma and 
tan (5 can still be used for parametrising the general Higgs sector in the R-parity violating 
MSSM. M\ is the mass of the lightest CP-odd Higgs boson in the R-parity conserving 
MSSM; as such, it is used here as a parameter. m\ is used to denote the physical tree- 
level mass of the lightest CP-odd Higgs in the R-parity violating MSSM (the convention 
adopted is that masses in the RPC case, parameters in this model, are denoted by M, and 
the masses in the RPV model are denoted by m). 



3.1 CP-even neutral scalar masses and couplings 

The Lagrangian after spontaneous gauge symmetry breaking contains the terms 



C D 



Re hi Re ulo Re i>Li ) 



/ Rehl \ 



2 

EVEN 



Rez/Lo 
V Re i>Lj J 

As such, the scalar CP-even Higgs squared mass matrix becomes 

[ cos^ /5Ml + sin^ /3M| -i sin2/3(Ml + M|) -Bj \ 
-isin2/3(M| + M|) sin^ /3M| + cos^ /5M| 5^tan/3 
-Bi Si tan/? M?S^ 



(3.2) 



EVEN 



V 



(3.3) 



where 



(Ml). + I 



- cos 2/5M| , 



(3.4) 



are in fact the sneutrino physical masses of the RPC case. It is important here to notice 
that the top-left 2x2 sub-matrix is identical to the RPC case, for which the Higgs masses 
are given by 



Kh = ^(Ml + Ml±^iMl + MIY - AMlMl cos^ 2/3 



(3.5) 



and will be used as parameters in the RPV model. 
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The matrix ()3.3|) always has one eigenvalue which is smaller than M|. This may be 
proved as follows: one first observes that the upper left 2x2 submatrix of ()3.3|) . call it A, 
has at least one eigenvalue smaller than or equal to Mf . Then using the Cour ant-Fischer 
theorem [21] of the linear matrix algebra one proves that, for one flavour, the eigenvalues 
of the 3x3 matrix A^Iven' ^'^^ interlaced with those of A. This means that the matrix 
Ai^vEN with i = 1 has at least one eigenvalue smaller or equal than M|. Repeating this 
procedure twice, proves our statement. Furthermore, it is interesting to notice that in the 
region where tan /? 1, the eigenvector (sin /?, cos /?, 0, 0, 0)""" corresponds to the eigenvalue 
with mass approximately M|. Notice that this is the same eigenvector as in the RFC 
case which corresponds to the Higgs boson which couples almost maximally to the Z-gauge 
boson. 

Lepton flavour violating processes have not been observed as yet and therefore, bearing 
in mind cancellations, the parameters Bi tan f3 have to be much smaller than min{M\^ ^f)- 
To get a rough estimate, consider the dominant contribution from neutral scalars and 
neutralinos in the loop [22], 

~ T~a ^3 ~ 1 , 3.6 

with rfi = max(MAi Mi) and B ~ OiBi). This shows that 

5itan/3 1.2- 10^3 



0.1% . (3.7) 



With this approximation, it is not hard to flnd a matrix Zr which rotates the fields into 
the mass basis, such that 

Zr'^TWIvenZr = diag m^o,m^o, (m?_^)i , (3.8) 
with m?o being the lightest neutral scalar mass and 



/ . cos(/3— «) cos aB, , sin(/3— a) sin a_B, \ 

/ cosa sma - .^.h^m-.mI) + co^m'^-M-J) \ 

cos(/3— o) sin oBj . sin(/3— q) cos oBj 

— Sm a COS a r-r.. R(Mi-Mi\ + r.r..RfMi-Mi \ 



I cospP^Bj cos/3Pf Bj ^ , 

\ cos(/3— a) sin(/3— a) *J / 

where there is no sum over i and (iWj, M|, Mfj) are defined in ()3.4I3.5|) . In addition. 



(3.9) 



Ml+Ml urr Ml cos^ 2/5 - M2 
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(the common convention is to choose < < 7r/2 and — 7r/2 < a < 0). The mass 
eigenstates of the RPV model are therefore given by 



R /,2 • R ~ , , COS PP!^B, \^ ^ 

COS a Re n.n — sm a tie ulo + 77, r Re vii 

cos(p — a J 

Q -r. , 2 -r. ~ I COS PP/^Bi _ 

H ^ sm a Re + cos a Re ulo + ( . ^ ^ 7 ) Re z/^j 



sin(/3 — a) 



cos{P — a) cos aBj sin(/5 — a)Bj sin aBj > 2 
' Re Air 



cos 



P{Mf - Ml) 



+ 



cosP{Mf-Mfj 



+ 



cos(/9 — a) sin aBj sin{(3 — a) cos aB 



cosp{Mf-Ml) 



+ 



cos P{M] - Mjj] 



Re ulo + Re i>Li 



with corresponding masses, 



(3.11) 



M, 



2 M| cos2 2/3 - Ml 



(M^^ - Ml) cos' P Ml - Ml 



5? 



B' 



M6 cos4/5' 



(3.12) 



, . M| cos^ 2/5 - 



Ml) cos'' M? 



B} 



i=l 



Ml 



+ 0{ 



B' 



M6 cos^f3' 



(3.13) 



[m 



cos^/i) 



Mf - M| cos2 2/5 



- Ml {Ml + M|) + Ml Ml cos' 2(3 



+ 0{- 



B' 



-,)■ 



(3.14) 



M6 cos'^p' 

The above expressions, are useful in relating the masses of the neutral scalars in the RPC 
and RPV case in the valid approximation i?tan/5 << min{Ml,Ml). They are presented 
here for the first time except the mass in ()3.14|) which agrees with Ref. [15]. We note 
here that these formulae are not valid if some of the diagonal entries in the mass matrix 
are closely degenerated - in such case even small Bi terms lead to the strong mixing of 
respective fields. However in many types of calculations (e.g. various loop calculations) 
one can still formally use such expansion - in the final result one often gets expressions of 
the type which have a well defined and correct limit also for degenerate masses, 

even if the expansion used in the intermediate steps was, in principle, wrong. 

It is interesting to note that the rotation matrix U defined in (12. (ij) , although explicitly 
calculated in this article, does not appear to all the neutral scalar vertices. For example. 
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the vertices of the CP-even neutral scalars with the gauge bosons read as^, 



>CvvH = I (cos (3Zn2s + sin (3Znis) 

2 cost^,,, 



+ -g^Mw (cos ^Zn2s + sin W^^ 



(3.15) 



where -ffg^i ... 5 are the Higgs boson fields, (i>+)i, (z/+)2, (t'+)3 respectively. From 

()3.9|1 and £vvh above, it is easy to see that the light Higgs boson coupling to the vector 
bosons (y = Z, W), is proportional to sin(/5 — a) as in the RPC case^. In fact, the coupling 
sum rule, 



E2 2 
9vVH0 — 9vv<f> ' 



(3.16) 



s=l 



valid in the RPC case for s = 1,2, persists also here, where gvvHO are the couplings 
appearing in ()3.15j) and gvv4> the corresponding coupling appearing in the Standard Model. 



3.2 CP-odd neutral scalar masses and couplings 

For the CP-odd case one finds, 



£ D — ^ Im hi Im Uiq Im Un ^ Ai 



2 

ODD 



/ Imhl \ 



-^ODD 



v 



J 



(3.17) 



where the CP-odd mass matrix reads, 

I cos2/?Ml + esin2/?M| i sin 2/3(Ml - ^M|) Bj \ 

isin2/?(Ml -^M|) sin^ /?Ml + ^ cos^ /3M| S^- tan/3 , (3.18) 

Bi Bi tan (3 

and ^ is the gauge fixing parameter in gauge. In fact, by using an orthogonal rotation 



V 



/ sm(3 — cos /5 \ 
cos sin P 

V 01/ 



(3.19) 



^Note that the matrix Z defined in 12.21|l has nothing to do with neither Zr nor Za defined in this 
section. 

^We foUow the conventions of Ref. [23]. 
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we can always project out the would-be Goldstone mode, of the CP-odd scalar matrix and 
thus 

/ iMl \ 



Ml 



cos P 



\ 







(3.20) 



cos /J I / 



Under the approximation of small bilinear RPV couplings [see Eq. fj3.7|) ]. a solution is 
determined for the matrix Za which rotates the fields into the mass basis, such that 



Za^A^oddZa = diag 



2 2 I 2 \ 



\ 



( sin [3 cos 13 

— cos [3 sin (3 
n 

^ cosl3{M'f~M\) 

with the mass eigenstates given by 

~ sin /? Im /ig — cos (3 Im z/^q i 

A° ~ cos /3 Im /iq -|- sin (3 Im u^q + 
-Bi ^ , n Bj tan /5 



Mf-Ml 
Bj tan /3 



(3.21) 



(3.22) 



5,; 



cos/3(M2-M2; 



-Im , 



M2-M2 



with corresponding masses. 



Ml 



1 ^ 

.«2 X] 



cos2 Mf-Ml 
B? 



+ 0{ 



B' 



+ 0{- 



M6 cos4/3' 
B^ . 



(3.23) 



(3.24) 



(3.25) 



(Ml - M2) cos2 /3 ' "'M6cos4/?^ 

The coupling of the Z gauge boson to the CP-even and CP-odd neutral scalar fields is given 
by 

3 



c 



ZHA 



-W2 

2cw 



{Pho —PAo)fi ZjH2+a)s ZA{2+a)p — ZruZaIp 



a=0 



Z^H^.Al , (3.26) 



where the four momenta p'^oyP^o incoming and the fields 5 correspond to 

A", (z>_)i, (z>_)2, (z>_)3 respectively. One may check that the coupling Z — — 
derived from ()3.2(i|l is proportional to sin(a — P) as it should be. 
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4 Positiveness and stability of the scalar potential 



4.1 Positiveness 



In general, one should inspect whether all physical masses in the CP-odd and CP-even 
sector are positive. For that, all diagonal square subdeterminants of mass matrices should 
be positive. One can easily check that both CP-odd and CP-even mass matrices in ()3.3I3.18|) 
respectively, lead, in the rotated basis, to the same set of conditions, 

> with. = 1,2,3 and Ml > -1- ^ . (4.1) 

i=l * 

Using the form of in ()3.H) . the last equation can be rewritten in the form 

i?o>tan/?X^^. (4.2) 

i=i « 

Excluding some very singular mass configurations, the above conditions are rather trivially 
fulfilled if one takes into account the bound of Eq. ()3.7|) . 



4.2 Stability 

The question of whether the potential is stable, i.e. bounded from below, is far more 
complicated. In most cases the quartic (/)— )term dominates and there is no problem. The 
only exception being when the fields follow the direction |/i2p = Ylt=o \ ^Li\^- In such a 
case, one should check whether the remaining part of the potential is positive along this 
direction. 



Denoting R = yYli=o\^i\'^ ^2 = -Re"**^, where is a free phase, and using 
Eqs. fj2. 2012. 2513. 4|) . one can write down the scalar potential along this direction in the 
vanishing snueutrino VEV basis as 



Bo 



Kcutral = J^^loho + [Mi + Bo cot (3] U^iU^i 

sin COS l) 

+ Bi tan/3 (P^o^^Li + ^^lo^^L) " ^« ('^La^o + H.c.) 
= uIQu-l- {B'^ hRe-'''' + R.c.) . (4.3) 
where the real symmetric matrix Q is 

= I ^"i ^ I . (4.4) 

' BjtanP [Mf + Bo cot P]6i, ' 
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Finding the stability conditions for the potential ()4.3|) is difficult, it depends on nine real 
variables (4 moduli and five phases of the fields). To simplify the problem, we perform one 
more field rotation to the basis in which the matrix Q is diagonal. This can be done, in 
general, by numerical routines (routines where already used in calculating the vanishing 
sneutrino VEV basis, and therefore, finding the stability conditions for the general scalar 
potential always has to involve some numerical analysis). We thus define the matrix P, 
Vl Pi>L, as 



ptQP = diag(Xo,Xi,X2,X3) . 



(4.5) 



In fact, Q is real, so we can choose P to be real orthogonal. Also, we denote -D^ = BaPafS- 
Obviously, the rotation P preserves the value of R = |/i2|. 

The potential becomes: 

3 



[Xa\ha\' - DaR (^L^e"^* + H.C.)] 



(4.6) 



a=0 



where Xq has to be positive, otherwise for = along the direction i/^j = Im z/^^q = the 
potential Koutrai = |Re z>^qP[Xo — i3)osign(Re 2^2.0)] ^^^^^ least for one direction along 

the Keuj^Q axis. In fact the condition on is X^ > 2\Da\- Thus our first conclusion is 
that the matrix Q has to be positively defined. One can write down appropriate conditions 
in the same manner as for the scalar mass matrices; comparing with Eq. ()4.H) . it can be 
observed that this condition is automatically fulfilled if relation 1)4.11) holds. 

With Xa positive, one can write down the potential as: 

2 3 ^2 



3 

a=0 



-.Re 



i(f) 



a=0 



(4.7) 



To further simplify the problem, denote z>^q, 
and (pa are free phases. Then 

3 



<t>oc) ^ where > are field moduli 



R' 



va=0 



'X. 



-Ur 



where R 



E 



i=0 



u. 



R 

^ Phases ( 



Da 




(4.8) 



can be adjusted independently of Ua- The 



worst case from the point of view of potential stability, the smallest first term inside the 



parenthesis, occurs for D^e 



i4>a 



\Da\. Denoting further = Ua/R, < < 1, one can 



reduce our initial problem to the question whether the function 



E 

a=0 



\Da\ 

JX. 



Ef^ 

a=0 



(4.9) 



0=0 
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depending now on four real positive parameters, is non-negative on the unit sphere 
Sa=o^a = 1. In general such problem can be solved numerically using the method of 
Lagrange multipliers. For Xi > Xq — Dq, the minimum occurs for 

where A can be found numerically as a root of the following equation: 

For smaller Xi, the minimum is realized for = for one or more values of i and requires 
analysis of various special cases. Having found the correct minimum, to prove the stability 
of the potential one needs to show that the function g at the minimum is non-negative. 

As shown in Eq. ()3.7|) . Bi terms and thus also Di terms are usually very small. In 
this case one can set approximate, sufficient conditions for the stability of the potential, 
without resorting to solving Eq. (jHHj), numerically. Denote D = Y^i=i^l ^^"^ Xmin = 
mm(Xi,X2,X3). Then, using the inequality < \/E?=i A^^/ELi^I = - 4^ 

one has 



9{ea) > Xoel + X.^U^-el) + {X,-X^,n)e!~2\Do\eo-2D^l-el. (4.12) 

Terms {Xi — Xmm)e? are always non-negative. The worst case being when the vector 
(^1; ^2, £3) is along the minimal Xi axis, where these terms vanish. Other terms are rotation 
invariant in the 3-dimensional space (ei, €2, 63), so Eq. ()4.12j) is equivalent to finding param- 
eters Xq, Xmin, -Do, D for which the expression ()4.13p . depending on just one real variable, 
is positive: 

g'ieo) = Xoel + X^,„(l - el) - 2\Do\eo - ^oJl-el > . (4.13) 



Analysis of ()4.13j) is further simplified by one more approximation, justified for small D: 

gXe^) > XQel + Xmin{l-el)-2\DQ\eQ-2D . (4.14) 

The rhs of Eq. ()4.14|) is now trivial. Following approximate conditions for the stability of 
the potential can be summarized as follows: 



Xmin range 


Stability requires 


Xmin > Xq — Dq 


Xq > 2\Dq\ + 2D 


< Xmin < Xq — Dq 


{Xq — Xmin){Xrain — 2i5) > D^ 



Both conditions are sufficient, but not minimal - we have made some approximations and 
there may be parameters which do not fall into either of the categories above, and yet still 
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give a stable potential. For example, if Xq = Xi = X2 = X3 = X, one can easily derive 
the exact necessary and sufficient condition for potential stability as X > 2\J + D"^, less 
strict than X > 2{\Do\ + \D\) which would be given by the table above. 
For complementary work the reader is referred to Ref. [24]. 



5 Conclusions 

In this letter we present a procedure for calculating the rotation matrix which brings the 
neutral scalar fields of the general R-parity violating MSSM onto the vanishing sneutrino 
VEV basis where they develop n zero VEVs, with n being the number of flavour generations. 
In doing so, we have made no assumption about the complexity of the parameters. We 
consider the case of n = 3 generations, but our approach immediately applies to other 
cases, apart from obvious modifications of the form of Z matrix defined in ()2.21I2.22|) . As 
a byproduct of basis change, we prove that the tree level MSSM potential does not exhibit 
any form of CP-violation, neither explicit nor spontaneous. Consequently, the neutral 
scalar fields can be divided into CP-even and CP-odd sectors with the 5x5 neutral scalar 
squared mass matrices, taking a very simple form with only small RPV masses sitting 
on their off diagonal elements. We can thus expand along small RPV masses and find 
analytic approximate formulae which relate the RPC and the RPV neutral scalar masses. 
Furthermore we also find, that in general there is always at least one neutral scalar field 
with mass lighter than Mz which couples maximally to the Z-gauge boson in the case of 
large tan j3 and large M^. Our procedure for finding the rotation matrix U has been coded^ 
and is numerically stable. 

In the end, we are aiming to construct the most general MSSM quantum field theory 
structure resorting neither to R-parity violation nor to other approximations. This will be 
useful for examining the phenomenology of the MSSM as a whole. The convenient choice 
of the basis for the neutral sector found in this paper is a ffist step towards this direction. 
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